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^—i ■ Abstract 
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I/"") ! The Belle Collaboration observed an enhancement called A (4630) in the A+A c mass dis- 

tribution using initial state radiation. We demonstrate that the enhancement could be consistent 
with the ip'J o(980) molecular picture of the F(4660) taking into account the A+A~ final state 
interaction. To test the hypothesis that the A(4630) and F(4660) are the same molecular state, 
we give predictions for its spin partner, the rj' c fo (980) molecule. High statistic measurements of 

^^ ■ the B decays into the KA^A~ and Krj' c n + 7r~ are strongly recommended. 
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The recently observed open and hidden charmed hadrons have stimulated many studies. They 
challenge our current knowledge of hadron spectroscopy, and provide us with an opportunity to un- 
derstand non-perturbative QCD better. Among these hadrons, the Y(4660) was observed by the Belle 
Collaboration in the ijj'tt + tt~ mass distribution using the technique of initial state radiation (ISR) (H- 
The mass and width were reported to be 4664 ± 11 ± 5 MeV and 48 ± 15 ± 3 MeV, respectively. This 
structure is very special because it was neither observed in e + e~ — > 7/s , i?vr + 7r~ J/i/j (H, nor in the 
mass distributions of a charmed and anti-charmed meson pair in the final states of electron-positron 
collisions (3HU. Furthermore, the 7r + 7r~ invariant mass spectrum shows a single peak at the high end, 
i.e. towards the mass region of the scalar meson /o(980). In Ref. HI it was argued that these facts 
may be naturally explained in terms of a hadronic molecular picture, i.e. by ip' fo(980) being bound 
together in an S- wave, while they would challenge other explanations (6j|71[8l. 

More recently, the Belle Collaboration reported another structure, called A (4630), in the A+A~ 
invariant mass distribution in e + e~ — > 7isrA+A~ J9). The reported mass is 4634^^g MeV, and the 
width is 92^24^21 MeV, consistent with the ones reported for the Y(4660) within two sigma. Based 
on the tetraquark picture, both structures were proposed to be of the same origin in Ref. 00, however, 
there is no general consensus on this issue yet (see e.g. the discussion in the short review |[i"0l ). In 
this paper, we shall show that they could also be understood as the same state within the ij/ fo(980) 
hadronic molecular picture, and discuss how this hypothesis can be tested in future experiments. 

In the -i///o(980) hadronic molecular picture, one may expect naively that the bound state would 
decay mainly through the decays of the unstable /o(980), and hence into the ip'tni, and the peak in 
the 7T7T invariant mass spectrum close to the /o(980) mass region appears naturally. While the latter 
statement is correct, the former one needs to be scrutinized. The mass of the Y(4660) is higher than 
open charmed and anti-charmed meson thresholds, and the A^A^ threshold. If the binding energy 
e = M^i + rn,/ (98o) — -^Y(4660) is ver Y small, the coupling of the bound state to its constituents 
determined by the equation JT2l 

|-=4(M^+m /o(980) ) 2 y^, (1) 

with /i the reduced mass of the ip' and /o(980), is small, and so is the partial width T(Y(4660) — >• 
i/j'tttt). On the other hand, the open charm channels have larger phase space, and might have larger 
partial decay widths. In fact, there is a well-known example — the /o(980) decays mainly into two 
pions which have plenty of phase space although it can be understood as a KK bound state lfi~3l[i"4l . 
In this paper, we shall assume that the A+A~ is the dominant open charm channel and study the 
implications of this assumption. This means we shall assume the total width of the Y (4660) is given 
by the sum of the partial widths into the tp'nir and A+A~, i.e. 

Ttot _ 3 [,/,'Tr+Tr-] „[A+A^] ,-. 

1 y — -1 y +ly , \4) 
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where the factor 3/2 in front of Ty is from isospin symmetry. 

The line shape of the Y(4660) is given by its spectral function 

py(m) = y y\ > — _ ; (3) 

M 2 -M Y + liy(M) 

convoluted with phase space, where My is the mass, ry*(M) is the energy-dependent total width, 
and Ily(M) = Ily(M) — Re [fly (My)] is defined as the self-energy with the real part subtracted at 



the mass |[T5l . The self-energy for arbitrary values of M is given by a dispersion integral (for further 
details, see Ref. HI) 

thr 

where M t h r denotes the relevant physical threshold. In Ref. (5], only the decays Y — )• ip'irir(KK) 
were considered. In order to check whether or not the structure observed in the A+A" mass distri- 
bution is consistent with the y(4660) observed in the il)'ir + it~ , one needs to include the contribution 
of the A+A~ in the total width Ty t - For that, a simple Lagrangian for the y(4660)A+ A~ coupling, 
which is assumed to be in an S wave, is used 

Cya c a c = -5yA c A c A c 7 M y M A c , (5) 

with gy a c a c a dimensionless coupling constant. Then the cross sections for e + e~ — > 7/5ij^ / 7r + 7r _ 
and e + e~ — > jisrA.^A~ are simply given by the corresponding parts of the spectral function of the 

y(4660) 

<r(A+A 
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where T Y ' n n andly c are the partial decay widths of the Y (4660) into the ^j'tt + 7t^ and A+A~ 
channels, respectively. The overall normalization constant N is the same for both processes since both 
structures were observed by the Belle Collaboration in the ISR processes. 

Since the Y (4660) has the quantum numbers J PC = 1 , it couples to the A+ A~ system in an S- 
wave, specifically to the 3 S\, and, therefore, the impact of the final state interaction (FSI) is expected 
to be large. In principle, the situation is comparable to J/ip decays with the proton-antiproton channel 
in the final state where FSI effects are known to play a rather important role |[i~6l [i~7l [i~8l [i~9l l20l I2T1 
22l . Unfortunately, there is no direct experimental information on the interaction between charmed 
and anti-charmed baryons. Thus, we have to resort to a model of the A+A~ interaction for taking 
into account FSI effects. Here we adopt the potential presented in Ref. |[23l . which was derived 
using SU(4) flavor-symmetry arguments, and compute the Jost function J{M) for this interaction. 
Multiplying the reaction amplitude with the inverse of the latter quantity, also known as enhancement 
factor, is practically equivalent to a treatment within a distorted- wave Born approximation Il24ll25l . 
The width of Y (4660) ->■ A+A~ is then given by 

yA e A c p A | . M l ' 

\J(M)\ 2 6tt I ^ M 2 



rl A ^](M) = ^4-^ (l + 2^ A 



Y ( M ) = i ,;\X C n £z I 1 + 2 ^7T 1 9 ( M ~ 2M A C ). (7) 



where M\ c is the mass of the A c , p = */M 2 /4 — M\ is its three-momentum in the rest frame of 

the F(4660), and 6 is the step function. In the calculations of Ref. 123], the function l/\J(M)\ 2 for 
the 3 £i channel decreases from about 2 at zero momentum to 0.3 at p ~ 500 MeV, and then slowly 
approaches unity only at very high momenta. In our calculations, we parameterize 1/\J(M) | 2 up to 
p ~ 500 MeV with the following function 

1 p 2 + b 2 



\J{M)\ 2 p 2 + cp + a 2 
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Figure 1: The FSI enhancement factor l/\J{M)\ 2 (dashed line) and the quantity p/\J(M)\ 2 (solid 
line) as a function of the excess energy e = M — 2M Ac . The latter curves are normalized arbitrarily, 
(a): the 3 S\ channel; (b): the 1 5q channel. 



with the parameter values being a = 247.7 MeV, b = 1390.4 MeV, c = 387.3 MeV, and d = 
0.0677. Then we set d = 1, which may always be done because such a normalization can be absorbed 
into a redefinition of the coupling constant gy a c a c , so that the remaining factor approaches unity 
asymptotically, and provides an enhancement to the amplitude close to the threshold. In Fig. Q] (a), 
the FSI enhancement factor in the 3 5i channel as well as this factor times the two-body phase space 
are shown as a function of the excess energy e = M — 2Ma c . Note that the central value of the peak 
observed by the Belle Collaboration in the A+A~ mass distribution is about 90 MeV above threshold, 
hence it cannot be due to the FSI enhancement solely, as may be seen from the figure. An opposite 
claim was made recently in Ref. ll26l . 

Using Eqs. ©, we perform a simultaneous fit to the cross sections of both processes. For simplic- 
ity, we assume that there is no background. Then there are three free parameters: the normalization 
constant N, the mass of the y(4660), My, and the Y(4660)A C A C coupling constant gy a c a c - The best 
fit gives 



N 



237J& 



My = 4662.5+^ MeV, g YAcAc = 0.7 ± 0.1, 



(9) 



with x 2 /d.o.f. = 1.4. The uncertainties quoted above are only from the fit, and do not include 
an estimate of the systematic uncertainty of the procedure. In doing the above fit, we chose to use 
M^i as given by the PDG ll27l and the central values of the parameters for the /q(980) measured 
recently by the KLOE Collaboration in the best fit Kl shown in Table 4 in Ref. ll28l . i.e. we used 
m,f = 976.8 MeV, gf K+K~ = 3.76 GeV and gj on + n - = —1.43 GeV The comparison of our best 
fit with the experimental data is presented in Fig. |2j cf. the solid lines. Also shown are the results for 
the case without the A C A C FSI (dashed lines), which were obtained with the same parameters except 
for the coupling constant. We use gyA c A c /\J{My)\ as the coupling constant for the case without FSI 
such that it coincides with the FSI modified coupling at the mass of the Y(4660). From the A+A~ 
mass distribution, one immediately sees the enhancement effect of the FSI on the cross section close 
to the threshold. From the best fit, we obtain the partial widths of the Y (4660) 



r(F(4660) -> t/Ztt+O = 8 MeV, T(y(4660) -> A+A~) = 93 MeV, 



(10) 
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2: The A+A~ and ip'n + 7r~ invariant mass spectra. The data are taken from the Belle measure- 
The solid curves are the results of the best fit, and the dashed curves are the results with FSI 
omitted. 



and their ratio is 



T(Y -► A+A-) 



11.5. 



(11) 



The ratio is smaller than the central value 24.8 extracted in Ref. (8j[TTl considering also an interference 
of the resonance with a polynomial background. In Ref. |[8l the authors also treated the X(4630) and 
the y(4660) as the same state, however, in this case as a compact tetraquark. 

At this stage, we want to emphasize that the FSI obtained from the model of Ref. (23\ is afflicted 
with sizeable uncertainties. However, it incorporates all essential features one expects from a realis- 
tic FSI, specifically it is generated by solving a scattering equation and it includes effects from the 
presence of annihilation channels. Therefore, it should be sufficient to give an illustration for the FSI 
effect in the problem at hand. The A+A^ interaction of Ref. 1231 contains two parts — an elastic 
part based on meson exchange and derived via SU(4) flavor symmetry, and an optical potential to 
simulate annihilation processes. In order to check in-how-far changes in the FSI influence our results 
we varied the strength of the optical potential by factors in the range from 1/2 to 2. It turned out that 
these variations only have a marginal effect on the resulting invariant mass distributions from the best 
fit. 

It should be clear that what we discussed above is only a possible scenario. The fact that one 
can obtain a combined fit of the A+A~ and the TTirtp' channels also in the molecular picture does not 
prove that the X(4630) and the Y(4660) are the same state. Observables should be found to further 
support or disprove this hypothesis. In this context, it is important to investigate the spin partner. 
Heavy quark spin symmetry in any case predicts the existence of a spin partner, however, the scenario 
outlined implies some very specific properties of that spin partner with respect to its mass and decay 
properties, as we will discuss now. 

In Ref. ll29l . based on heavy quark spin symmetry, we predicted the presence of an 77^/0 (980) 
bound state, called Y v , as the spin multiplet partner of the ip'fo bound state. The mass of the Y v 
should satisfy 

M Yv = M Y(4mo) - (Afy - Mrf e ) (12) 

to a high precision. Using the best fit value for the Y(4660) mass given above and M„/ = 3637 ± 
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Figure 3: Predictions of the Y T) line shapes in the r]'n + ir and A+ A c in arbitrary units. The solid and 
dashed curves represent results with and without FSI, respectively. 



4 MeV E71 . one gets My v = 4613 ± 4 MeV where the uncertainty is dominated by the one from 
the r]' c mass. Based on the same formalism as above, the line shape of the Y v in the 7]' c tt + tt~ and the 
A+A~ may be predicted. Heavy quark spin symmetry indicates that the coupling of the Y n to the 
A+A~ has the form, cf. Eq. ©, 



£y„A c A c = Wya c aAc7 5 Y v A c 



(13) 



with the same coupling constant as the Y(4660). 

In Fig. |3l the predictions for the Y n line shapes in the rj'ir + TT~ and A+ A^T channels are shown in 
arbitrary units, however, with the relative normalization fixed. With the FSI, now in the 1 Sq partial 
wave and calculated again from the A+A~ model of Ref. ll23l . shown in Fig.[T](b), the predicted line 
shapes are given by the solid curves, while the ones without FSI are given by the dashed curves. The 
Y n mass is only about 40 MeV higher than the A+ A~ threshold, as a result the width of the Y v is much 
smaller than that of the Y(4660), and thus the line shapes are much narrower. The partial widths for 
decay into the r]' c ir + ir~ and the A+A~ channels are 8 MeV and 22 MeV, respectively. The ratio 



T(Y V -> A+A-) 

T(Y V -)■ T//7T+7T-) 



2.7 



(14) 



is much smaller than the one for the Y(4660) as a result of smaller phase spaces. Furthermore, the 
effect of the FSI is not so significant anymore. We expect that within other models for the spin partner 
of the Y(4660) the discussed properties, especially the mass and the ratio of Eq. (fl"4l) . will be very 
different. 

In summary, taking into account the A+A~ FSI, we found that the X(4630) may be described as 
the same state as the Y(4660) in the t///o(980) bound state picture. One notices that there should be 
other open charm decay channels, such as decays into charmed and anti-charmed mesons. We checked 
that an additional constant width from other possible decay channels of less than 30 MeV may still 
be accommodated. In principle, a polynomial background as in Ref. (H allows one to improve the 
fit. Also possible interferences of the A(4630) or y(4660) with other resonances, such as highly 
excited ip resonances, could have an impact on the analysis. However, neither of these effects is under 
control quantitatively given the current quality of the experimental data. Hence, in our analysis, we 



refrain from considering them to reduce the number of parameters. Within the molecular picture for 
the Y(4660), the presence of a Y v with a mass given by Eq. (fT2l as the spin partner of the ip'fo(980) 
bound state is almost unavoidable, since the spin-dependent interactions are highly suppressed by 
1/m^, with m c the charm quark mass J29ll . Other models of the Y(4660) should also provide a spin 
partner, but most probably with a different mass and different decay patterns. Thus, in order to test 
the molecular picture it is important to search for the Y n experimentally, for instance in the decays 
B^ — > t] i c K ± it + it~ which is expected to have a large branching fraction |29l . 

At last, we want to mention that a related observation was made by the BaBar Collaboration in the 
reaction B~ -)■ A+A~K ~ OJl. They observed a structure at 2931 ± 3 ± 5 MeV in the A+A~ mass 
distribution. In the paper, the A+A~ mass distribution is also provided, where one can see clearly two 
peaks. The measured branching ratio of the decay B~ — > A+A~ K~ is of order 10~ 3 11301 , which is 
several orders higher than the naive expectation 10~ 8 since this three-body decay is color-suppressed 
and with a small phase space 0T1 . In Ref. Oil Cheng et al. showed that the high suppression could 
be diminished, if there was a narrow hidden charm state with a mass of order 4.6 — 4.7 GeV or a 
charmed baryon, which was assumed to have J p = l/2 + , coupled to the A+A~. We notice that the 
positions of the double peaks coincide with the masses of the Y(4660) and the predicted Y v . However, 
they could also be due to a charmed baryon H c with J p = 3/2 + — we found that a J p = l/2 + H c 
baryon, as used in Ref. OH . cannot describe the double peak structure in the A+A~ mass distribution. 
Also some interference of a charmed baryon with the charmonia is possible. Better data with higher 
statistics, especially better Dalitz plots, would be very helpful in illuminating the situation. 
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